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ABSTRACT 

We investigate the influence of photoevaporation of protoplanetary discs on the 
final distribution of exoplanets semi-major axis distances. We model giant planet mi¬ 
gration in viscous discs affected by photoevaporation driven by either pure EUV or soft 
X-ray radiation (XEUV). We show that the final exoplanet distributions are strongly 
dependant on the choice of the photoevaporation model. In particular, we find that 
XEUV is more efficient than pure EUV radiation at parking planets at approximately 
1-2 AU distance from their central star, hence roughly reproducing the observed peak 
in the exoplanets semi-major axis distributions. We note however that a more quan¬ 
titative comparison with the observations is hindered by the oversimplified treatment 
of planetary accretion, which severely affects migration rates. For this reason, caution 
should be used when using these models to constrain details of disc clearing and/or 
migration from the observations. Nevertheless our results indicate that disc disper¬ 
sal by photoevaporation may be the main driver of the features in the exoplanets 
semi-major axis distribution observed by recent surveys. 
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1 INTRODUCTION 

Statistical studies of the properties of extrasolar planets 
have recently become possible thanks to the discovery of 
over 1500 confirmed planets and almost 5000 total plan¬ 
ets, if the Kepler candidates are included (exoplanets.org, 
Han et al. 2014). In particular the study of the distribution 
of exoplanets semi-major axis could provide useful insights 
into the planetary migration and protoplanetary disc evolu¬ 
tion process (e.g. Alexander & Pascucci, 2012, from now on 
AP12). The occurrence of planetary migration is reasonably 
well accepted for high mass planets, which migrate inwards 
due to tidal interactions with the protoplanetary disc (e.g. 
Goldreich & Tremaine 1980; Lin & Papaloizou 1986). Indeed 
a number of numerical models have been computed, which 
roughly agree with the observed distributions (e.g. Armitage 
2007; Alexander & Armitage 2009; AP12). The migration of 
planets larger than 0,5 Mj up follows the so-called Type II 
migration, which sets in when a gap in the protoplanetary 
disc is dynamically cleared by the planet itself. At this point 
the planet and the gap keep moving inwards due to the ma¬ 
terial which keeps flowing into the gap on timescales dictated 
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by viscous accretion of the protoplanetary disc. While gas 
is present in the disc, this process is unstoppable and irre¬ 
versible, meaning that without a disc clearing mechanism 
all giant planets would be expected to fall into the star at 
some point. Indeed, planetary population synthesis models 
(Benz et al. 2014) include simplified prescriptions for even¬ 
tually dissipating the disc. Photoevaporation of the gas disc 
by high energy radiation from the central star (e.g. Alexan¬ 
der et al. 2006; Gorti et al. 2009; Ercolano et al. 2009; Owen 
et al. 2010) provides a mechanism to remove the gas from 
the disc and cease migration. 

API2 construct numerical models of planet migration 
in evolving protoplanetary discs in order to investigate the 
effect of disc clearing on the final semi-major axis distribu¬ 
tion of exoplanets. They find that assuming a pure EUV- 
photoevaporation profile (Alexander et al. 2006) and a flat 
mass distribution of extra-solar planet the observed peak at 
semi-major axis distances of approximately 1-2 AU could 
be roughly reproduced (see Hagesawa & Pudritz for an al¬ 
ternative explanation for the origin of planetary system ar¬ 
chitectures). However, depending on the nature of the pho¬ 
toevaporating radiation (EUV, X-ray or FUV) the shape 
of the photoevaporation profile, as well as the photoevap¬ 
oration rate itself, change dramatically. For example, com- 
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bined X-ray and EUV (XEUV) photoevaporation rates are 
roughly 100 times larger than pure EUV photoevaporation 
rates (Owen et al. 2010), and XEUV photoevaporation rates 
as a function of disc radius are not so strongly concentrated 
at the gravitational radius as in the pure EUV case (see e.g. 
Armitage 2011). One could expect that these large differ¬ 
ences should result also in different final distributions of ex¬ 
oplanets semi-major axes, meaning that observed distribu¬ 
tions may provide a constraint on the dominant disc disper¬ 
sal radiation. Furthermore the initial planet mass function 
may also have an influence on the final exoplanet semi-major 
axes distribution. 

In this paper we perform numerical modelling to explore 
the influence of different photoevaporation models and of the 
initial planet mass function on the final distributions of ex¬ 
oplanets semi-major axis, in order to explore the possibility 
to use the latter to constrain the initial conditions and the 
nature of the photoevaporating radiation. In Section 2 we 
describe our numerical approach and present a number of 
tests of the new code we developed to this aim. In Section 3 
we present and discuss our results, while a summary of our 
findings is given in Section 4. 


2 MODELLING APPROACH 


Our modelling approach is similar to that described in AP12, 
based on the numerical model of Alexander & Armitage 
(2009). We summarise here the basic concepts for conve¬ 
nience. Each model includes a protoplanetary disc which 
evolves under the influence of viscosity in the disc and pho¬ 
toevaporation from its central star. After a given timescale 
tform, a giant planet of a given mass m p is assumed to form 
within the disc and undergoes Type II migration due to the 
torques from the gaseous disc. Angular momentum conser¬ 
vation dictates that the disc evolution is in turn affected by 
the planet. We do not model the formation of the planet, 
and assume that this is formed istantaneously in the disc. 
The disc-planet system evolves according to the following 
equation: 
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where the first term describes the viscous evolution of the 
disc (Lynden-Bell & Pringle 1974), the second term the mi¬ 
gration of the planet due to the torques from the disc (e.g. 
Lin & Papaloizou 1986) and the last term the mass loss due 
to photoevaporation (e.g. Clarke et al. 2001). Here E is the 
surface density of the disc, R the distance from the star 
in the plane of the disc, v the kinematical viscosity of the 
disc, M* the mass of the star and E w the photoevaporation 
profile. 

The term A is the rate of specific angular momentum 
transfer from the planet to the disc. We use the following 
form for A, which is a slight modification by Armitage et al. 
(2002) of the original one used by Lin & Papaloizou (1986): 


A (R, a) = 



if R < a 
if R > a 


( 2 ) 


where q is the mass ratio between the planet and the star, 
a is the semi-major axis of the planet orbit (assumed to be 
circular) and A p is given by 

A p = max(77, |R — a|), (3) 


where H is the disc scale height. 

Due to the disc torques, the planet migrates at a rate 
given by: 
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This approach for describing migration is computation¬ 
ally more efficient than more complicated treatments, yet it 
was shown to give comparable results (Takeuchi et al., 1996, 
Trilling 1998). In Type II migration, the planet behaves es¬ 
sentially like a test particle entrained in the viscous flow, and 
this is the fundamental reason why a more detailed descrip¬ 
tion does not yield significantly different results. Note, how¬ 
ever, that doubts have been cast recently on the correctness 
of Type II migration in hydrodynamical simulations (Duf- 
fell et al. 2014, Diiermann & Kley 2015). We ignored these 
recent results as they are still somewhat uncertain, but note 
that they might point to important effects to be included in 
future works. 

To solve equation [TJ we use standard numerical tech¬ 
niques. We discretise the equation on a grid of points eq- 
uispaced in R 1 / 2 . The viscous term is treated as a diffusive 
term (Pringle & Verbunt 1986). The planet torque is treated 
as an advection term, and computed after the diffusive term 
using operator splitting. We use the van Leer (1977) method 
to reconstruct the surface density at the cell boundaries. The 
photoevaporation term is integrated by simply removing a 
fixed amount of mass from each cell at every timestep. To 
prevent numerical problems, we use a floor surface density 
of 10 ~ s g/cm 2 . Following Armitage et al. (2002) we limit 
the maximum planet torque close to the planet for com¬ 
putational reasons; this has no consequences on the orbital 
migration rate. We use a resolution of 4000 points and our 
grid extends from 0.04 to 10 4 AU. This is higher than what 
commonly employes in similar treatments. 


2.1 Planetary accretion 


It is well known that gaps cleared by giant planets are not 
able to stop completely the mass flow from the outer disc 
(e.g., Lubow & D’angelo 2006). For this reason we follow 
partially Alexander & Armitage (2009) and include a pre¬ 
scription, calibrated on published results from hydrodynam¬ 
ical simulations, to include this mass flow. We use the same 
prescription from AP12, which was proposed by Veras & 
Armitage (2004): 
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( 5 ) 


Here, the efficiency e(M p ) is defined as the ratio between 
the accretion rate onto the planet and the accretion rate in 
a steady disc in absence of the planet. It follows that 


Mp = e(M p )M disc 


(6) 


Following AP12, the mass accretion rate through the gap is 
computed via 
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Mi nn£r = -- M p . (7) 

1 + e 

Note, however, that the sum of the two rates is not 
Mdisc, and therefore the correct rate to apply at the outer 
edge of the gap to enforce mass conservation is the sum of 
the two rates. We note that such parameterisation may be 
avoided by using the results of Lubow & D’angelo (2006), 
who showed via analytical arguments that Minner = M p /e. 
We chose however to follow AP12 for the sake of comparison, 
and note that given the current uncertainties in these rates, 
both prescriptions may carry large uncertainties. 

To apply this mass leakage through the gap computa¬ 
tionally, we set the cells outside the orbit of the planet to the 
floor value, until we have removed the mass that is flowing in 
the timestep. We then use the prescribed rate M p to increase 
the mass of the planet and we “unload” the mass coming 
from the rate Minner in the first cell inside the orbit of the 
planet. This algorithm follows from a private communica¬ 
tion with R. Alexander and was chosen for consistency with 
their results. Possible improvements would consist in remov¬ 
ing and adding the mass from a larger region, to smooth the 
effect of this procedure (e.g., Owen 2014). 


2.2 photoevaporation modelling 


Several photoevaporation models are present in the litera¬ 
ture, considering different irradiating fields. The EUV field 
was the first to be explored (Hollenbach et al. 1994) through 
analytical calculations. The original estimates were then re¬ 
fined through hydrodynamical simulations by Font et al. 
(2004). Finally, Alexander et al. 2006 have shown the im¬ 
portance of including the direct illumination from the cen¬ 
tral star once a hole has been opened in the disc. We follow 
AP12 and assume a constant EUV flux $ = 10 42 photons 
per second in our models. We note however that this pa¬ 
rameter is currently poorly constrained, and recent obser¬ 
vational results (Pascucci et al. 2014) imply that this may 
have been significantly overestimated. The total mass-loss 
rate with this flux is 4 x 10~ 10 Mq yr~ 4 (Alexander et al. 
2006). The mass-loss profile E(r) can be found in the ap¬ 
pendix of Alexander & Armitage (2007). We distinguish be¬ 
tween a diffuse profile, which is the relevant one for most of 
the time while there is a full disc, and a direct profile, which 
is the relevant one after a hole has opened in the disc and the 
inner edge of the outer disc is directly exposed to the EUV 
radiation. A criterion is needed to switch between the two. 
AP12 apply the switch when the surface density in the inner 
disc falls below a critical value E c = 10~ 5 g cm -2 . Follow¬ 
ing a private communication with R. Alexander, we monitor 
the surface density at two different locations in the disc, at 
0.1 AU and 1 AU from the star, and switch between the two 
profiles when both have fallen below the critical value. The 
exact criterion assumed is important, since when the direct 
profile is applied the disc is dispersed in a short time-scale, 
and the planet migration is essentially stalled. 

In addition to the original EUV model used by AP12, 
we employ here also the XEUV model of Owen et al. (2010, 
2011, 2012). In this case the mass-loss rate is given by the 
following relation (Owen et al. 2012): 


M = 6.25 x 10 -9 Mq yr _1 


Lx 


1.14 


(8) 


where Lx is the X-ray luminosity of the star. The X-ray lu¬ 
minosity function is well characterized by observations (e.g. 
Preibisch et al. 2005 and Guedel et al. 2007), and the me¬ 
dian is 1.1 x 10 3 ° erg s~ 4 for solar type stars. Note that the 
jmass-loss rate corresponding to this X-ray luminosity is al¬ 
most 100 times higher than the EUV one. The shape of the 
E(r) function is parameterised in appendix A of Owen et al. 
( 2012 ). 

Finally, a photoevaporation model based on FUV radi¬ 
ation has been presented by Gorti, Hollenbach & Dullemond 
(2009). However, this work does not include a hydrodynam¬ 
ical solution for the wind, meaning that the resulting rates 
are at this time still uncertain and will not be considered 
here. 


2.3 Population synthesis 

We construct a population synthesis of disc-planet systems 
by randomly sampling a number of initial conditions, which 
are summarised for the two models in Table [T] The two dif¬ 
ferent photoevaporation models that we consider have very 
different mass-loss rates. This would lead to differences in 
the lifetimes of the disc if the same initial conditions are 
assumed. To reconcile the population synthesis with the 
observed disc lifetimes, the two different photoevaporation 
models need thus to assume different initial disc properties. 
The EUV model, which has a lower photoevaporation rate, 
assumes a slighly less massive, faster evolving disc than the 
XEUV one. Note that, since both the disc viscosity and the 
disc masses are poorly constrained, both models are a rea¬ 
sonable description within the current uncertainties. In all 
of our models we assume the initial surface density profile 
to be 

E(R,t = 0) = E 0 (^)exp(-A) ; (9) 

where Eo is a normalization constant that we fix by imposing 
a given total mass of the disc, and Ro is the initial size of 
the disc. We assume that the viscosity profile follows v oc R. 
This is convenient from the theoretical point of view, as it 
corresponds to a mildly flared disc and a constant value of 
a throughout the disc. In what follows we will speak of the 
normalization of the viscosity profile through the viscous 
time t v {R) = R 2 /3u(R) oc R. 

In the EUV model, we follow AP12 in varying the ini¬ 
tial disc mass. We assume a log-normal distribution, with a 
median of lCU 1 ' 5 Mq and a 3er dispersion of 0.5 dex. The vis¬ 
cous time is 5 x 10 4 yr at 10 AU. Note that in Id models there 
is a degeneracy between a, the Shakura-Sunyaev parameter, 
and H/R, the aspect ratio of the disc. The viscous time is 
the relevant quantiy that controls the evolution of the disc, 
and not the a parameter. For reference, we note that assum¬ 
ing a value of a = 0.01 and an aspect ratio of H/R = 0.0333 
at 1 AU gives the assumed viscous time. We follow AP12 in 
assigning the properties of the formed giant planets. Rather 
than modelling the formation process, we prescribe statisti¬ 
cally the properties of the formed planets, as the purpose of 
this exercise is to assess the implications of disc dispersal, 
rather than building a fully consistent planetary population 
synthesis model. We assume that all planets form at 5 AU, 
with a mass assigned from a uniform distribution between 
0.5 and 5 Mj up . We will later correct statistically this as- 
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Figure 1. Comparison Spock against Ivanov analytical solution 
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Figure 2. Comparison Spock against AP12 


sumption. The other important variable is the formation 
time, which we also draw randomly from an uniform distri¬ 
bution between 0.25 Myr and t c = l/3t v (3Md/{2t^M w )) 2 ^ 3 
(Clarke et al. 2001, Ruden 2004), which is the time of disc 
clearing. 

In the XEUV model, we follow Owen et al. (2011) in 
assuming that all discs start with the same initial condi¬ 
tions, and the only difference in their evolution is due to 
the spread in X-ray luminosities. Here we assume an ini¬ 
tial mass of 0.07 Mq, and a viscous time of 7 x 10 5 yr at 
18 AU. We use the X-ray luminosity function from Guedel 
et al. (2007), obtained for the Taurus cluster, to prescribe 
the X-ray luminosity of the star, which sets the mass-loss 
rate. The properties of the formed planets are assigned as 
before. 

2.4 Code tests 

We have tested independently the different components of 
our code, Spock. The viscous evolution part is the easiest 
to check, since analytical solutions (Lynden-Bell & Pringle 
1974) are readily available. No significant deviation has been 


observed. We tested the implementation of the photoevapo¬ 
ration profiles together with the viscous evolution by check¬ 
ing that the lifetimes of the discs (without planets) reported 
by Alexander et al. (2006,2009) and Owen et al. (2010) are 
correctly recovered. 

The planet migration part can be tested against an ex¬ 
isting analytical solution (Ivanov 1999). A comparison is 
shown in Figure [T| It should be noted that the analytical 
solution employs a different treatment for the planet, and 
therefore an exact overlap between the two is not to be ex¬ 
pected. Our result is similar to other codes (e.g., see Figure 
3 in Lodato & Clarke 2004 for a similar test). 

Figure[2]shows a comparison run of our code with AP12 
for a disc of mass 10~ 1,5 Mq. The two codes show reasonable 
agreement, particularly when considering that a number of 
relatively unconstrained parameters have to be introduced 
in the model, which may influence the final result. We have 
however no way of telling if those differences may become 
more serious for planets of different masses, as API2 runs 
at different masses are not available. We fear that this may 
be indeed the case when we compare the final distributions 
of exoplanets we obtain from our code for the same initial 
conditions parameter range as AP12, as shown by the red 
dotted lines in the left and right panels of Figure [3] Both 
for the high (right) and low (left) mass exoplanet distribu¬ 
tions, we do not find the strong peaks reported by AP12 
at semimajor axis values between 1 and 2 AU, indeed our 
distributions for the EUV case are relatively flat. Without 
a dedicated benchmarking effort of the two codes, the ex¬ 
act cause of this discrepancy is unknown, but it probably 
lies in the (often arbitrary) choices of a number of param¬ 
eters describing the redistribution of mass due to accretion 
on the planet and the leakage through the gap, which in 
turn has influence on the migration rate of the planet it¬ 
self. The other very sensitive part of the code is the switch 
between the direct and diffuse profile, which rapidly dissi¬ 
pates the disc and is therefore very important to treat cor¬ 
rectly. Until progresses are made in the physical description 
of planetary accretion and migration, the predicting power 
of models like those presented here or in AP12 is severely 
reduced and hence their use as diagnostics of disc clearing, 
planet migration and planetary accretion should be treated 
with caution. 

Nevertheless the main aim of this work is not affected 
by these uncertainties, as we wish to investigate the influ¬ 
ence of the photoevaporation model and the exoplanet mass 
function on the final exoplanet distribution. These aims can 
be achieved regardless of the uncertainties inherent to the 
modelling that are highlighted above, as long as selfconsis¬ 
tence within the same code is observed. 


3 RESULTS 

Following AP12, in our further discussion we divide our 
planetary sample into two mass bins: the low mass planets 
comprising planets of masses from 0.5 Mq to 2.5 Mq and 
a high mass bin comprising planets with Mq > 2.5 Mq. We 
explore in what follows the importance of the chosen photo¬ 
evaporation model and initial planet mass function assumed 
on the final semi-major axis distribution of exoplanets. 
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Model 

EUV 

XEUV 

Disc mass 

median 10 —1,5 Mq, log-normal with variance 0.5 dex 

0.07 M 0 

Scaling radius 

10 AU 

18 AU 

Star mass 

1 Mq 

0.7 Mq 

Viscous time at 10 AU 

10 4 yr 

4 x 10 5 yr 

X-ray luminosity 

/ 

X-ray LF from Taurus, Guedel et al. (2007) 

EUV flux 

10 42 phot/yr 

/ 


Table 1. Disc properties in the two models 


3.1 Effects of photoevaporation models on the 
final semi-major axis distribution of 
exoplanets 

The left and right panels of Figure [3] show a comparison 
of the final semi-major axis distribution of the low and high 
mass exoplanets formed in discs under the influence of either 
EUV or XEUV photoevaporation. The distributions shown 
in the Figure were obtained assuming a flat initial mass func¬ 
tion of exoplanets. As mentioned already in the previous sec¬ 
tion, the distributions obtained for discs under the influence 
of pure EUV photoevaporation are rather flat. The XEUV 
case shows on the contrary peaked distributions, with the 
peak being at semi-major axes values roughly comparable 
to those observations. It is worth keeping in mind, when 
comparing the XEUV to the EUV case that the accretion 
properties of the discs in these two cases are very different. 
In order to obtain similar disc lifetimes, the XEUV case re¬ 
quires a larger accretion rate to disc mass ratio than the 
EUV case, this leads to faster accretion of mass onto the 
planets in the XEUV case, and the overpopulation (about 
800 planets) of the high mass bin (right panel in Figure [3]) 
compared to the EUV case (about 150 planets). Once a 
planet reaches a given mass, the flow of material through 
the gap is strongly reduced, aiding the onset of a PIPE- 
like dispersal (Rosotti et al. 2013), i.e. inner disc dispersal 
by photoevaporation, which stops further migration of these 
planets. The final result is that a much larger number of high 
mass planets are retained in the XEUV case, compared to 
the EUV. Indeed both models start with 250 planets in this 
mass bin with the EUV case ending with only 150 planets 
in this mass bin, compared to the 800 planets in the XEUV 
case. The peak seen in the XEUV case is indeed due to this 
parking effect brought by the fast inner disc removal. 

It is clear from the above that the final distribution of 
extrasolar planets semi-major axis is seriously affected by 
the choice of photoevaporation model. As discussed above 
the main driver of the strong differences, however, is the 
different total mass loss rates predicted by the two models, 
which lead to differen accretion properties that have to be 
assumed in the two cases in order to preserve the correct 
disc lifetimes. The question then remains as to whether the 
photoevaporation profile itself has an influence on the final 
distribution of exoplanets semi-major axes. To answer this 
question we have constructed a further population synthesis 
where we keep all mass accretion rate parameters the same 
as in the EUV case, but apply the XEUV mass loss profile, 
normalised to a total mass loss rate of 4 x 10 -10 Mq yr -1 , 
which corresponds to the pure EUV model. This ensures 
that the correct disc lifetimes are preserved while isolating 
the effect of the mass-loss profile on the final distributions. 


The results of this academic experiment are shown in 
Figure [4] For low mass planets very little difference is seen 
between the EUV and XEUV profile cases, however the sit¬ 
uation is dramatically different in the high mass bin. Here 
many more planets are parked at large radii when an XEUV 
profile is used instead of an EUV for the same disc condi¬ 
tions. The disc lifetimes are very similar in the two cases, 
since the total mass loss rates have been set to the same 
(EUV) value, in fact the onset of dispersal by photoevapo¬ 
ration is even delayed in the XEUV case, given that for the 
same mass loss rate the profile is broader. The broader pro¬ 
file for the XEUV case leads also to migration rates for high 
mass planet in the XEUV profile case that are slower than 
in the EUV case. This happens because the XEUV is effi¬ 
cient at removing mass from the disc out to several tens of 
AU, while EUV-driven photoevaporation is mostly concen¬ 
trated at the gravitational radius (see e.g. Armitage, 2011). 
Therefore a planet that migrates inwards in an XEUV pho- 
toevaporated disc will be subjected to weaker gas torques as 
the surface densities in the disc will be lower over a broader 
range of disc radii. So by the time the inner disc is removed 
by photoevaporation, planets in an XEUV photoevaporated 
disc will have migrated less than planets in an EUV photo- 
evaporated disc. Indeed those in an EUV photoevaporated 
disc are free to fall inwards undisturbed at least until they 
reach the gravitational radius of the EUV irradiated disc 
(~1AU which is indeed where we see a small pile up in the 
EUV case). It is also worth keeping in mind at this point 
that the migration rate through the disc is only slowed down 
if E * af < M p , otherwise migration proceeds at the viscous 
time irrespective of the disc mass. This explains why we only 
observe this effect in the high mass bin. 


3.2 Effects of the initial mass function on the final 
semi-major axis distribution of exoplanets 

The exact shape of the initial exoplanet mass function is 
still somewhat uncertain. In the above tests and in AP12 
a flat mass distribution of planets was assumed. Given the 
different migration rates and the different influence of higher 
mass planets on the transport of the material across the gap, 
which determines the support of the inner disc, one could 
suspect that the shape of the initial mass function of exo¬ 
planets should also influence the final exoplanets semi-major 
axis distribution. We have compared the final distributions 
obtained for the XEUV and the EUV cases assuming a flat 
initial planet mass distribution and a simple 1/M distribu¬ 
tion and find the final distributions to be qualitatively simi¬ 
lar in the two cases. This is shown in Figure 5 for the XEUV 
case. The pure EUV case yielded similar results for the two 
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loglO(Final planet location / AU) loglO(Final planet location / AU) 


Figure 3. Histogram of exoplanet semi-major axes distribution for low - 0.5 Mq < M p < 2.5 Mg- (left panel) and high mass stars - 
M p > 2.5 Mq- (right panel). The black and red dotted lines are for the XEUV and EUV cases, respectively. A flat initial mass distribution 
of planets is assumed. 
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Figure 4. Histogram of exoplanet semi-major axes distribution for low — 0.5 Mq < M p < 2.5 Mq- (left panel) and high mass stars - 
M p > 2.5Mq- (right panel). The black and red dotted lines are for the XEUV and EUV profile cases, respectively. In this case only the 
photoevaporation profile changes (see text for details). A flat initial mass distribution of planets is assumed. 


initial mass distributions and it is not shown here. Both for 
the EUV and XEUV case the assumption of a 1/M initial 
planet mass distribution yield to a slight overpopulation of 
planets parked at slightly larger radii. The plots were con¬ 
structed considering the mass of the planets at the end of 
the simulation, comparing the mass of the planets injected 
at the beginning of the simulation increases slightly the im¬ 
portance of the bins at larger semi-major axis distances, but 
the results remain qualitatively the same. 


4 CONCLUSIONS 

We have used a numerical approach to investigate the in¬ 
fluence of internal photoevaporation of protoplanetary discs 
on the final exoplanet semi-major axis distance. We found 
the final distributions of exoplanets semi-major axes to be 
extremely sensitive to the choice of photoevaporation model 


in the gaseous disc. This implies that exoplanets semimajor 
axis measurements, which are becoming statistically signif¬ 
icant with the more recent surveys, may be theoretically 
used to constrain the details of disc dispersal. Our results 
show that an XEUV photoevaporation profile naturally re¬ 
produces the observed peak of planets residing at roughly 
1-2 AU distance from the central star. The XEUV photo¬ 
evaporation model is more efficient at parking planets at 
larger distances from the central star than a pure EUV pho¬ 
toevaporation models. This happens thanks to the ability of 
the XEUV model at removing gas from an extended region 
in the disc in contrast to the pure EUV model, in which the 
mass loss rate is concentrated at the gravitational radius. 
An XEUV photoevaporated disc experiences a lowering of 
the gas surface density at a range of disc radii, which slows 
down migration particularly for the larger planets. We cau¬ 
tion, however, that while our study aims mainly at theoret¬ 
ically comparing the influence of difference photoevapora- 
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Figure 5. Histogram of exoplanet semi-major axes distribution for low — 0.5 Mq < M v < 2.5 Mg- (left panel) and high mass stars 
— M p > 2.5 Mq- (right panel). The black and red dotted lines are for the XEUV case comparing, respectively, a 1/M initial mass 
distribution of planets with a flat initial mass distribution. 


tion models on the final semi-major axis distribution of ex¬ 
oplanets, a more detailed comparison with the observations 
remains difficult, due to a number of poorly constrained pa¬ 
rameters on which the theoretical models rely. In particu¬ 
lar, migration rates, which are crucial to the problem, are 
severely affected by details of the planet accretion formalism 
employed. In any case, our results also highlight the need of 
employing more precise treatments of photoevaporation in 
planet population synthesis models. 
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